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Abstract 

In a recent report, a robotic spacecraft mission is proposed for the purpose of collecting a small asteroid, or 
a small part of a large one, and transporting it to an orbit in the Earth-Moon system. Such an undertaking 
will require solutions to many of the engineering problems associated with deflection of an asteroid that poses 
a danger to Earth. In both cases, it may be necessary for a spacecraft to approach an asteroid from a nearby 
position, hover for some amount of time, move with the same angular velocity as the asteroid, descend, 
perhaps ascend, and finally arrest the angular velocity of the asteroid. Dynamics and control in each of 
these activities is analyzed in order to determine the velocity increments and control torque that must be 
provided by a reaction control system, and the mass of the propellant that will be consumed. Two attitude 
control algorithms are developed, one to deal with synchronizing the spacecraft’s angular velocity with that 
of the asteroid, and the other to arrest the asteroid’s angular velocity. A novel approach is proposed for 
saving fuel in the latter case. 
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1. Introduction 

The National Aeronautics and Space Administration (NASA) is currently investigating a conceptual 
robotic mission to collect a small Near-Earth Asteroid (NEA) or a small boulder resting on the surface of a 
large NEA, and transport it to an orbit in the Earth-Moon system. An object with a diameter of roughly 
7 m is of interest. Such a mission would lead to greater scientific understanding of NEAs and serve as 
preparation for planetary defense. This paper is concerned with dynamics and control of the spacecraft 
trajectory and attitude in the terminal phases of a rendezvous with the object, as well as the attitude of the 
spacecraft after its payload has been secured. 

The objective of NASA’s investigation is to consider in further detail the mission proposed in Ref. [1], 
Figure 1 illustrates aspects of the mission that are especially relevant to a study of dynamics and control. 
An asteroid, A, is to be collected by a spacecraft, B, using a container that is essentially axisymmetric. In 
Fig. 1, b 3 is a unit vector fixed in B and parallel to the container’s axis of symmetry. A few assumptions 
that strongly influence the direction of study can be made based on the configuration of the container. 
First of all, placement of the container’s opening requires B to advance in the direction of b 3 as it collects 
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A. Secondly, as long as the diameter of the container exceeds the largest dimension of A, collection can 
be accomplished without restricting the relative orientation about t >3 between A and B. We assume that 
in order to prevent damage to the container, two requirements must be imposed: as B approaches A, b 3 
must remain collinear with a line fixed in A, and there must be no relative angular velocity between A and 
B when the two make contact. Finally, we assume that the angular velocity of B relative to an inertial 
reference frame must be made to vanish after A is collected. This is necessary, for example, to orient solar 
arrays properly or to aim thrust in the direction needed to return to the Earth-Moon system. 

Small asteroids, which are more numerous than large ones and collide with Earth more frequently, could 
with sufficient warning time be deflected by using such a spacecraft. The spacecraft could also be used to 
collect a small asteroid for the purpose of deflecting a large one, as proposed in Ref. [2] . 

Throughout what follows, the spacecraft B is regarded as a rigid body. As discussed in Ref. [1], the mass 
of B is 18,000 kg. The bus (without the container) is axisymmetric, with a diameter of 2.7 m and a height 
of 5.9 m. For the purpose of this study, we presume the container is light in comparison to the bus, and 
mass is distributed uniformly in the cylindrical bus. In that case I r , central principal moments of inertia of 
B (that is, principal moments of inertia relative to the mass center of B) for unit vectors b r (r = 1, 2, 3), are 
determined to be I\ = / 2 = 60416 kg-m 2 , and / 3 = 16403 kg-m 2 . Four pods, each containing four Reaction 
Control System (RCS) thrusters, are mounted to the bus on struts as shown in Fig. 5 of Ref. [1] (p. 23). 
Each pod has a radial distance of 2 m from the cylinder’s axis of symmetry, and is 0.3 m from the cylinder’s 
end. Each thruster can apply 200 N of force, and has a specific impulse I sp = 287 s. The position vector 
r B Pr from the mass center of the bus, B*, to each of the four pods P r ( r = 1,2, 3, 4), and directions of 
thrust, are reported in Table 1. 


Table 1: RCS Thruster Configuration 


Pod 

Position (m) 

Thrust Directions 

Pi 

r B * Pl = 2b ! - 2.65b 3 

±b 2 , ±b 3 

Pi 

t b*p 2 = _ 2 b 1 _ 2.65b 3 

±b 2 , ±t>3 

P 3 

r s * p 3 = 2b 2 - 2.65b 3 

±t>i, ±b 3 

Pi 

v B*p i = _2b 2 - 2.65b 3 

±bi, ±b 3 


Whereas the foregoing description of B can be given quantitatively, a corresponding concrete description 
of A is not available because the particular object to be collected has not yet been identified. Moreover, 
there is a paucity of accurate information about 10-m asteroids, and such information is not easily obtained. 
Consequently, to make progress in analysis we must fabricate our own asteroids and assign each of them 
a heliocentric orbit, physical dimensions, mass, distribution of mass, and rotational motion. In practice, it 
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will be important to design the spacecraft to be able to deal with a healthy amount of uncertainty in these 
and other asteroid parameters. 

We consider the process of collecting A to consist of the following steps. B leads or trails A in its orbit at 
a nearby position that is stationary relative to A; B then makes an initial approach and begins to hover with 
respect to A. Hovering continues for an appropriate amount of time. The orientation of B is changed so 
that the opening of the container is facing A, and angular velocity of B relative to A is eliminated. Descent 
takes place as b 3 remains collinear with a line fixed in A. The object of interest is collected in the container 
and made firmly attached to B; if A is a large asteroid, B must ascend from its surface. The final step 
involves arresting the angular velocity of the rigid body formed by B and its payload, denoted simply by B 
for convenience. Trajectory and attitude dynamics and control involved in each of these steps are examined 
in detail in Secs. 2-6 of the paper. Concluding remarks appear in Sec. 7. 

2. Initial Approach 

Our analysis of dynamics and control begins with the initial approach of the spacecraft from a point 
relatively close to the asteroid. In the interest of simplicity, the heliocentric orbit of the asteroid, A , is taken 
to be a circle of radius R, as shown in Fig. 2. Throughout the paper, A is treated as a rigid body, and 
the gravitational force exerted by A on the spacecraft is neglected. Rotational motion of A is modeled as 
torque-free, before it is captured. 

Analysis is facilitated by introducing a right-handed set of mutually perpendicular unit vectors Iq, L 2 , 
and L 3 , fixed in a local-vertical-local-horizontal reference frame L, and directed as shown in Fig. 2. L 3 is 
in the direction of the position vector from A * , the mass center of A, to S 7 the Sun. Li is in the direction 
of the asteroid’s orbital velocity or, more precisely, the velocity of A* in a Newtonian reference frame N. 
Together, Li and L 3 define the plane of the heliocentric orbit. 



Figure 2: Asteroid in Circular Heliocentric Orbit 
The angular velocity N u> L of L in N is given by 

N uj L = -n L 2 ( 1 ) 

where H is a positive constant, and L 2 is normal to the asteroid’s orbit plane. The most general rotational 
motion of A that can be considered is tumbling; we expect, however, that it will prove difficult to carry 
out and present results of studies of approach, hover, descent, and ascent associated with a tumbling body. 
Consequently, we conduct a preliminary investigation by taking the angular velocity l uj of A in L to be 

L u > A = uj Li (2) 
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where w is a positive constant. As is well known, the proposed rotational motion is not stable unless 
N u ) A = u;Lii — HLo is parallel to the central principal axis of inertia of A associated with the largest central 
principal moment of inertia, or the smallest. In what follows, we concern ourselves with intervals of time 
during which l uj a satisfies Eq. (2). Furthermore, the equatorial plane of A is regarded as the L 2 -L 3 plane, 
the direction of the North pole is L 1; and the direction of the South pole is — Li- 

Before it commences to approach A, the spacecraft is taken to be stationary in L at a point P, trailing 
the asteroid over its South pole as shown in Figs. 2 and 3. The approach trajectory ends in conditions that 
permit hovering to begin; in other words, the velocity of the spacecraft is such that it is fixed in reference 
frame (rigid body) A. In Fig. 3, Q belongs to a set of points on the surface of a sphere, D, whose radius 
exceeds the largest dimension of A. We investigate trajectories to points on the surface of D. 



Figure 3: Hovering Stations and Collision Avoidance Constraint 

The total velocity increment, AU, needed for each trajectory is minimized with an interior-point opti- 
mization algorithm that allows the trajectory to be constrained in order to avoid a collision with A. The 
constraint can be described as a simple geometric shape that envelops A; in this study, the shape is taken to 
be a sphere, denoted by C in Fig. 3. The algorithm is applied to translation governed by linearized equations 
of motion relative to a circular orbit, sometimes referred to as the Euler-Hill equations (see Ref. [3], p. 111). 
Here, we neglect terms in those equations representing perturbing forces acting on the spacecraft, such as 
gravitational attraction of the asteroid and solar radiation pressure. 

The results that follow are obtained with R = 0.9 astronomical units (AU), which corresponds to f2 = 

l. 3 x 10 -5 deg/s. Each trajectory departs from P, 50 m behind A*. Time of flight is arbitrarily specified 
as 2600 s (43.3 minutes), which we believe to be a reasonable duration. The points on the surface of D at 
which each trajectory ends are at various latitudes and longitudes with respect to A; the radius of D is 20 

m. The radius of C, the constraint or “keep-out” sphere, is set to 10 m. w takes on values of 0.5, 1, 5, 10, 
and 60 rev/hr (rph). 

Under the aforementioned conditions, trajectories ending over the southern hemisphere of A require 
two impulses: one to depart P, and one to yield a velocity such that the spacecraft can begin to hover 
at the desired point fixed in A. The radius of C is small enough that some trajectories ending over the 
northern hemisphere can be accomplished with two impulses. Other trajectories to stations over the northern 
hemisphere require three impulses. 

As one might expect, for stations at a given latitude the required AU is virtually independent of the 
longitude. The results, therefore, can be presented for various values of latitude and w, as seen in Fig. 4. 
Latitudes of —90°, 0°, and 90° are associated, respectively, with points over the South pole (— Li), equator, 
and North pole (Li) (see Fig. 2). The left side of Fig. 4 contains plots of AU; the plots on the right side 
show the propellant mass, Am, obtained from AU via the rocket equation Am = mo[l — e~ Av ^ a Isp ' > ] with 
spacecraft mass m 0 = 18, 000 kg, J sp = 287 s, and g = 9.8 m/s 2 . 


4 



-90-75-60-45-30-15 0 15 30 45 60 75 90 



-90 -75 -60 -45-30-15 0 15 30 45 60 75 90 




latitude (deg) latitude (deg) 


Figure 4: Ay and Propellant Mass for Approach Trajectories 


Results for stations directly over the poles are independent of w; the two-impulse trajectory to the South 
pole requires AV = 0.023 m/s and Am = 0.15 kg, whereas the three-impulse trip to the North pole requires 
AE = 0.075 m/s and Am = 0.48 kg. For low spin rates, u> = 0.5 and 1 rph, the trajectory to the North pole 
is more expensive than one to a point over the equator. On the other hand, reaching the hovering condition 
over the equator becomes the most expensive maneuver at the higher spin rates. 

In the case of a small asteroid for which L u> satisfies Eq. (2), it is only necessary to consider approaching 
a point on the South or North pole. It is, however, worthwhile to consider approaching points at various 
latitudes because a suitable boulder may be resting on the surface of a large asteroid at any location. 


3. Hover 

The approach trajectories studied in Sec. 2 put the spacecraft in position to hover; however, propulsive 
force is required to maintain hover. A hovering spacecraft traces out in L a circle whose center is on the 
line that passes through A* and is parallel to (see Fig. 2). The plane of the circle is parallel to the L 2 -L 3 
plane. Because the stations of interest lie on the surface of D centered at A* (see Fig. 3), the radius p of 
the circle can be written as 

p= H cos A (3) 

where H is the radius of D , and where A is the latitude of the station. Thus, for a spacecraft regarded as a 
particle Q hovering at a station that is initially in the plane of the asteroid’s orbit and on the same side of 
A * as the Sun, the position vector from A* to Q can be expressed as a function of time t , 

r A Q = H sin ALi + H cos A(— sin LutL 2 + cos LutL 3 ) (4) 

From this relationship one can obtain quantities appearing in the Euler-Hill equations, 

x = r _i ^ • L 3 = 77 sin A y = r A< ^-L 2 = —77 cos A sin tut z = r A Q • L 3 = H cos A cos tot (5) 
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Figure 5: A V and Propellant Mass for Hovering 


and the Euler-Hill equations can then be written as 

x — 2ilz = 2HQu>cos Asinwf = f • Li (6) 

y + f l 2 y = H(uj 2 — fl 2 ) cos A sin Lot = f • L 2 (7) 

z + 2ttx — 3 Vl 2 z = —H(lu 2 + 3fl 2 ) cos A cos tot = f • L 3 (8) 

where f is the resultant force per unit mass acting on Q, other than the gravitational attraction of the Sun. 
The smallest value of co presently under consideration is 5 x 10 -2 deg/s, whereas f 1 = 1.3 x 10 -5 deg/s. 
Therefore, in comparison to terms involving to 2 , one can neglect terms involving and f l 2 . Thus, after 
multiplying Eqs. (7) and (8) by sinwf and — cos wf, respectively, and adding the results, one obtains 

Hu; 2 cos X = f • (sinwtL 2 — coswfLa) = — f • a 3 (9) 

where unit vector a 3 is fixed in A and has the same direction as L 3 at t = 0. This result tells us that 
hovering can essentially be accomplished by applying propulsive force per unit mass f directed toward the 
spin axis of A , along a line passing through Q and perpendicular to the spin axis. It is therefore advisable 
for the orientation of the spacecraft to be such that a pair (or two) of thrusters can apply thrust in the 
direction of — a 3 . 

Letting a 3 and a 2 be unit vectors fixed in A such that a 3 = Li, and a 2 has the same direction as L 2 at 
t = 0, one can express the velocity increments needed to maintain hovering for one revolution of A as 

r‘t—2n/uj 

A V r = | f • k r \dt (r= 1,2,3) (10) 

Jt = o 

8 Oy r 

AEi = 811.77 cos A AV^ = — fl 2 H cos A AV^ = — (w 2 + fl 2 )H cos A 

to Ul 
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( 11 ) 


The total velocity increment AV = AVi + AV 2 + AV 3 (AV 3 w AV) required to hover for one revolution of 
A is plotted on the left side of Fig. 5 for various values of and w and latitude, A. The plots on the right side 
show the propellant mass, Am, obtained from AV via the rocket equation with spacecraft mass mo = 18, 000 
kg, I s p = 287 s, and g = 9.8 m/s 2 . Although Eqs. (11) are developed for stations that at t = 0 lie in the 
orbit plane of A , they are nevertheless applicable to all stations that share a common longitude because the 
integration in Eqs. (10) is performed over a complete revolution of A. Figure 5 indicates that pole-sitting 
can be accomplished without any effort whatsoever, which is to be expected because the polar axis of A is 
assumed to remain parallel to Li, and the stations over the poles are in the orbit plane of A at the same 
heliocentric distance as A*. Figure 5 also shows that hovering comes at a cost that increases proportionally 
with u). A comparison of Figs. 4 and 5 shows that hovering for one revolution is more expensive than the 
initial approach, unless the station is very near a pole. It is important to keep in mind that for, say, u> = 1 
rev/hr the associated curve for AV in Fig. 5 shows what is required to hover for one hour, whereas the curve 
for AV corresponding to u> = 60 rev/hr must be scaled by a factor of 60 to determine what is required for 
one hour of hovering. Consequently, hovering can require prohibitive amounts of propellant when to is large. 

4. Removal of Relative Angular Velocity 

Following a period of time during which hovering is performed, as discussed in Sec. 3, the spacecraft’s 
attitude and angular velocity are altered to prepare for collection of the object of interest. We assume 
the method of collection requires that there be no relative angular velocity between the asteroid and the 
spacecraft when the two make contact. Letting B denote the spacecraft, one can express this condition 
as A u> B = 0, where A u> B is the angular velocity of B in A. The condition can also be expressed as 
a uj b = n uj b — N u ) A = 0, where n lj a and n uj b are, respectively, the angular velocities of A and B in 
a Newtonian reference frame N . Clearly, a control system cannot be used to make N uj b equal to N u> ‘ 4 
unless N u> A is a known function of time. 

The reader will recall that we treat rotational motion of A in A as torque-free. When the orientation of 
A in N, and the values of N ui A • a r (r = 1, 2, 3), are specified at a particular time to, one can then produce 
a time history of N u> A in a number of ways. (In practice the orientation and N laJ a at f 0 are determined 
from measurements made with sensors placed on B.) The first of these is numerical integration of Euler’s 
dynamical equations together with a set of kinematical differential equations governing parameters suitable 
for describing the orientation of A in N. That is the approach we take in this work. Alternatively, closed- 
form solutions for n lj a • a r (r = 1,2,3) are available in terms of Jacobian elliptic functions when A is 
unsymmetric; that is, when it has three distinct central principal moments of inertia (see Ref. [4], Sec. 3.4). 
In the special case of two identical central principal moments of inertia, A is axisymmetric and its angular 
velocity can be described simply in terms of spin about the axis of symmetry, together with precession about 
the angular momentum vector (see Ref. [4], Sec. 3.1). 

As mentioned previously, we regard spacecraft B as an axisymmetric rigid body. We assume that 
operation of the capture mechanism requires the axis of symmetry of B to be collinear with a line of 
descent fixed in asteroid A , but that effective use of the mechanism does not require any particular relative 
orientation, about the symmetry axis, between B and A. 

In Ref. [5], Tsiotras and Longuski address the problem of controlling B while it performs simple rotational 
motion in N about b. 3 , a unit vector parallel to the axis of symmetry. The objective is for b 3 to have the 
same direction as h, a unit vector fixed in N. They use two parameters to represent the direction of n in B , 
and develop a globally asymptotically stable control algorithm. We have extended their work to develop a 
“spin-up” feedback controller that drives S 3 to the same direction as a, a unit vector fixed in A and parallel 
to the line of descent, and also drives a uj b to 0 . Whereas control is accomplished in Ref. [5] with torque 
applied in two directions, our objectives in general require three-axis control. 

A specific procedure is now described. We begin by identifying a line of descent fixed in A, and letting 
a be a unit vector parallel to this line. Predictions are made as to the direction of a in A at some future 
time t* , as well as the value of N u> 4 • a at t*. An initial reorientation of B is performed to bring S 3 to the 
predicted direction of a in A. Provided n uj b = 0 before the reorientation begins, the change in the attitude 
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Figure 6: Relative Angular Velocity and Alignment of Axis of Symmetry During Spin-up 


of B can be brought about slowly, by briefly firing RCS jets at the beginning and end of the maneuver, so 
that only a small amount of propellant is expended. (When the cost of hovering is high, it may be advisable 
to perform this slow maneuver before hovering begins.) Next, n uj B • b.3 is made to have the predicted value 
of N u> A • a. Finally, at time t = t* the feedback controller is employed to bring A u: B to 0 , and to give 
b 3 the same direction as a. Even though the motion of A is torque-free, control torque must be applied 
continuously to B to keep a uj b = 0 . 

Use of the feedback controller is demonstrated in an example involving an unsymmetric asteroid, taken 
to be an ellipsoid with semidiameters a = 5 m, b = 3 m, and c = 4 m, a uniform density of 2000 kg/m 3 , 
and therefore a mass of 502.65 x 10 3 kg. If we take ai, a 2 , and a 3 to be parallel to central principal axes 
of inertia of A, then the associated central principal moments of inertia are, respectively, J\ = 2.52 x 10 6 
kg-m 2 , J 2 = 4.13 x 10 6 kg-m 2 , and J 3 = 3.42 x 10 6 kg-m 2 . The angular velocity of A in N is characterized at 
t = t* = 0 by the values n oj a • a 3 = 0.6 deg/s, N uj a • a 2 = 0.6 deg/s, and N uj a • a 3 = 6.0 deg/s (60 rph). 
Unit vectors hi, n 2 , and n 3 fixed in N can always be chosen such that n r = a r (r = 1, 2, 3) at t = t* = 0. 
The line of descent is chosen to be parallel to a 3 . 

According to the procedure previously described, an initial reorientation of B brings b , 3 to the direction 
a 3 is predicted to have at t = t* = 0, and the required propellant mass is assumed to be negligible. At 
the conclusion of the reorientation, n lo b = 0 . Next, four RCS thrusters (see Table 1) are used to apply a 
couple whose torque is given by T = 1600 b 3 Nm; N u: B ■ b 3 reaches the required value of 6.0 deg/s in 1.07 
s, and 0.31 kg of propellant is consumed in the process. To emphasize that the relative orientation about 
b 3 between B and A can be arbitrary, we assume the angle between ai and b 3 (which is the same as the 
angle between a 2 and b 2 ) is 60° at t = t* = 0. At this time the feedback controller is employed, and the 
resulting motion of B relative to A is shown in Fig. 6. Plots of a uj b ■ b r (r = 1,2,3) appear respectively 
in the upper left, upper right, and lower left of Fig. 6; in the first two cases, the controller drives the errors 
to zero in approximately 200 s, whereas in the third case the controller is able to maintain an error of zero. 
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Figure 7: Control Torque and Propellant Mass for Spin-up 


The angle between t>3 and £3, which can be expressed as cos _1 (a3 • S3), is plotted in the lower right of Fig. 
6. The controller keeps this angle less than 4 deg, and drives the angle to zero within 200 s. 

Control torque must be applied continuously to B to keep a uj b = 0 for t > 200 s, as can be seen in 

Fig. 7. Time histories of r r = T • b r (r = 1,2,3), where T is the control torque exerted by the RCS on 
B, are shown respectively in the upper left, upper right, and lower left of Fig. 7. The torque needed to 
synchronize rotational motion of B with A is large because the asteroid is tumbling at high angular speed. 
The propellant mass A m consumed by the RCS thrusters, shown in the lower right of Fig. 7, reaches 31 kg 
after 1000 s. It is important to note that this expenditure of propellant gives t>3 the same direction as E13, 
which is parallel to the line of descent; additional propellant is required, however, to make B translate such 
that b3 and the line of descent fixed in A remain collinear. 

5. Descent and Ascent 

After a uj B has been brought to 0 and 63 has been given the direction of a line of descent fixed in A, as 
discussed in Sec. 4, B can begin to descend in the direction of b 3 towards A. When A is a large asteroid and 
a small boulder is to be collected from its surface, B must ascend after the payload is secured. Beginning 
from a condition of hover, descent is initiated with an impulsive velocity increment. Likewise, descent is 
terminated and B is brought to rest in A by an impulsive change in velocity. During the descent, thrust 
must be applied continuously to keep B from moving relative to A in any direction normal to the line of 
descent. Ascent is regarded as descent in reverse. 

In what follows we consider an asteroid moving as described in Sec. 2, and report AV and corresponding 
propellant mass Am for descent and ascent. These quantities are functions of latitude of the line of descent, 
time of flight, and lo [see Eq. (2)]. 

In examining descent, we regard A as sphere with a radius of 3.5 m. Each descent begins from a point 
20 m distant from A*, and ends at a point on the surface of A, 3.5 m from A*. In each case, the line of 
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descent is taken to be normal to the surface of A. The left side of Fig. 8 contains plots of AV; the plots 
on the right side show the propellant mass, Am, obtained from AV via the rocket equation with spacecraft 
mass mo = 18, 000 kg, I sp = 287 s, and g = 9.8 m/s 2 . Results correspond tow = 0.5 rph on the top, and to 
to = 60 rph on the bottom. (Although we also obtained results for to = 1, 5, and 10 rph, they are omitted 
in the interest of brevity.) The three curves in each plot correspond to values of latitude A = 0° (solid), 
±45° (dash-dot), and ±90° (dash); results are independent of the sign of A, as is the case with results for 
hovering shown in Fig. 5. 

Time of flight (the duration of descent) ranges from 100 s to 2700 s in the case of to = 0.5 rph, and 
from 5 s to 100 s in the case of w = 60 rph. The ranges are selected so as to illustrate the presence of 
minima in some of the curves; descending too quickly or too slowly can be expensive. Optimal times of 
flight for to = 60 rph are markedly shorter than those for to = 0.5 rph; although a rapid descent can save a 
significant amount of propellant, it may not be operationally prudent. For a given value of to , descending to 
the equator is more expensive than descending to any other latitude. Conversely, descending to either pole 
is the least expensive; in fact, descent along the polar axis is independent of to. For this reason, the curves 
associated with A = ±90° on the top and bottom of Fig. 8 are two parts of the same curve, which decreases 
monotonically. 

To study ascent, we treat A as a sphere having a radius of 250 m. Each ascent begins from a point on 
the surface of A, proceeds along a line of ascent fixed in A and normal to the surface, and ends at distance 
of 300 m from A* (in other words, at an altitude of 50 m). 

The left side of Fig. 9 contains plots of AV. Because larger asteroids are thought to rotate more slowly 
in general than smaller ones, we allow to to take on the values of 0.1 rph (top), and 1 rph (bottom). For 
convenience, the boulder to be collected is taken to be the object studied in Sec. 4; its mass is 502.65 x 10 3 
kg. Thus, the spacecraft and payload together have a mass of 520.65 x 10 3 kg, and this is the value of mo 
used in the rocket equation to obtain the curves plotted on the right hand side of Fig. 9. 

6. Removal of Inertial Angular Velocity 

After a small asteroid has been secured in the container of B, or after a small boulder has been collected 
and B has ascended from the surface of a large asteroid, n lo b is made to vanish. This is considered 
necessary, for example, to orient solar arrays properly or to aim thrust in the direction needed to return 
to the Earth-Moon system. We assume that the spacecraft and its payload now form one rigid body. For 
convenience B will denote the newly formed rigid body, henceforth regarded as unsymmetric in order to 
maintain generality. 

In this section we develop an asymptotically stable control algorithm for applying torque to drive n lo b 
to 0. No requirements are placed on the orientation of B in N at the time N to B reaches 0; the required 
orientation can subsequently be obtained by briefly firing RCS thrusters at the beginning and end of a slow 
reorientation that consumes a negligible amount of propellant. The controller is exercised with an example 
taken from Ref. [1]. 

Euler’s dynamical equations of motion for a rigid body B can be written as 

I • N a B = T - n lo b x I- N to B (12) 

where I is the central inertia dyadic of B , N a B is the angular acceleration of B in a Newtonian reference 
frame N , and T is the moment about B *, the mass center of B , of external forces acting on B. The 
rotational kinetic energy K of B in N can be written as 

I< = \ n lo B n lo b (13) 

Hence, the time rate of change of K is given by 

I<= n lo b I - N a B = N oo B • (T - n lo b x I • N uo B ) = N to B ■ T (14) 

Because the vector product of N to B and the vector I • N to B is necessarily perpendicular to N to B , the scalar 
product of N oo b and n lo b x I • N to B necessarily vanishes. 
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In order to drive n uj B to 0, we seek a strategy for applying control torque T that will dissipate rotational 
kinetic energy; that is, we wish to have K < 0. At the same time, we seek to apply torque sparingly so as to 
avoid excessive expenditure of propellant. These considerations lead us to an objective function J defined 
as 

J = u t t + | t t Rt (15) 

where w is a 3 x 1 column matrix whose elements are n uj B ■ 8 r , r is a 3 x 1 column matrix whose elements 
are T • b r (r = 1,2,3), R is a 3 x 3 positive definite matrix, and where a superscript T denotes matrix trans- 
position. Forming the matrix product uj t t yields the same result as forming the scalar product N u> B • T, 
namely the time rate of change of rotational kinetic energy. The matrix product |r T Rt is referred to as 
control energy. The matrix r for which J is a minimum can be found by setting the partial derivative of J 
with respect to r equal to zero, 

u t + t t R = 0 (16) 

whereupon we find 

r = - R (17) 

and 

K = -uj t R (18) 

Because R is positive definite (therefore, so is l? -1 ), AT < 0. Moreover, K = 0 if and only if u> = 0. Thus, 
according to Lyapunov theory (Ref. [6]), r as given by Eqs. (17) provides asymptotically stable control of 
the system governed by Eqs. (12). Numerical values of the elements of R can be adjusted to trade the use of 
control energy for dissipation of rotational kinetic energy: higher values result in smaller control torque and 
slower dissipation of energy, whereas lower values lead to larger control torque and faster energy dissipation. 

An example in Ref. [1] (p. 35) involves an asteroid modeled as a cylinder with a radius of 3 m, a height 
of 12 m, and a mass of 1.1 x 10 6 kg. The cylinder’s axis of symmetry is taken to be parallel to unit vector 
bi fixed in the spacecraft; thus, central principal moments of inertia of the asteroid for bi, b2, and b3 are, 
respectively, J± = 4.95 x 10 6 kg-m 2 , and J2 = J3 = 15.675 x 10 6 kg-m 2 . Keeping in mind that the height of 
the spacecraft is 5.9 m and the radius of the asteroid is 3 m, the position vector from the spacecraft mass 
center to the asteroid mass center is given by 5.95 b 3 m, and the central principal moments of inertia of B 
for bi, £>2, and £>3 are I\\ = 0.5673 x 10 7 kg-m 2 , K 2 = 1.6398 x 10' kg-m 2 , and A3 = 1.5691 x 10 7 kg-m 2 . 
Evidently B is unsymmetric as the values of Ai, A'2, and A3 are distinct. Angular velocity is taken to be 
n lj b = 6b 3 deg/s. It naturally occurs to an analyst to use the thrusters described in Table 1 to apply force 
in the form of a couple whose torque has the same direction as — b3, in which case the moment arms are 2 
m. The authors of Ref. [1] take such an approach and determine that bringing N uj B to 0 consumes 306 kg 
of propellant. (The value is 292 kg according to our calculations). 

We now show it is possible to use substantially less than 300 kg of propellant to arrest the rotational 
motion of B. After referring to Table 1 one can observe that, without the asteroid, r B Pr • 83 = —2.65 m. 
More importantly, with the asteroid in place, the value becomes r B Pr • 83 = —8.33 m. This moment arm 
is more than 4 times longer than 2 m; there is a corresponding improvement in efficiency if thrusters can be 
employed so as to take advantage of the longer moment arm. To bring the longer moment arm into play, one 
must first forgo application of thruster forces in couples. As will be demonstrated presently, in the current 
example one can ignore the consequences of applying thruster forces with a non-zero resultant. Moments 
about B * of thruster forces applied with the 8.33-meter moment arm are necessarily perpendicular to b3, 
whereas moments parallel to b3 are necessarily applied with 2-meter moment arms. Consequently, it might 
at first appear that the long moment arm cannot be used to advantage in this case because the central 
angular momentum of B in N has the same direction as 83. However, 83 is parallel to the axis associated 
with the intermediate central principal moment of inertia, A 3 , and rotational motion about such an axis is 
known to be unstable. It is possible, therefore, to exert torque briefly in a direction perpendicular to b3 
and intentionally induce instability to obtain tumbling motion; during the ensuing torque-free motion of the 
unsymmetric rigid body, the direction of angular momentum in B changes dramatically. When the angular 
momentum becomes nearly perpendicular to b3, the long moment arms can be brought to bear. (It is worth 
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keeping in mind that this approach could not have been used if, in this example, N uj b were parallel either 
to bi or to b 2 , because in each case the rotational motion would be stable. Furthermore, this approach 
cannot be used when B is axisymmetric because in that case a small change in angular momentum cannot 
lead to tumbling.) 

Figures 10 and 11 display the results of inducing tumbling, waiting for the central angular momentum of 
B to become perpendicular to b 3 , and then engaging the asymptotically stable control given by Eqs. (17). 
At time t = 0, the central angular momentum of B in N is given by H = I • n lo b = 1.6432 x 10 6 b 3 Nms. 
Torque given by T = 800 bi + 800 b 2 Nm is applied for 1 s. During the subsequent torque-free motion, 
n lo b • b 3 becomes zero three times as can be seen in the lower left of Fig. 10, with the last of these instances 
occurring at t = 1256.5 s at which time the controller is activated. Vertical dashed lines indicate t = 1256.5 
s in Fig. 10, and this instant is immediately obvious in Fig. 11. Control could have commenced the first or 
second time N u> B • b 3 reached zero; the third instance was chosen in order to illustrate periodic torque-free 
rotational motion that takes place after perturbing torque has been applied. The matrix R is taken to be 
diagonal with all three elements equal to 100. Rotational motion of B is essentially arrested by t = 2500 s, 
which is 1250 s after control begins. Time histories of T r = T • b r (r = 1, 2, 3) are shown respectively in the 
upper left, upper right, and lower left of Fig. 11. Maximum magnitudes of T\ and r 2 are larger than that of 
t 3 by about a factor of 10. Propellant mass Am consumed by the thrusters is shown in the lower right of 
Fig. 11 and it amounts to 108 kg, which represents a savings of 184 kg. 

The reduction in RCS propellant consumption may come at the expense of a disruption to the orbital 
trajectory because the resultant thrust force is not zero. This possibility, however, is in general mitigated 
to a large extent because the trajectory is normally controlled by means of a solar electric propulsion (SEP) 
system having an J sp = 3000 s, which greatly exceeds the 287-s J sp of the RCS. Another mitigating factor 
is that rotational motion may allow RCS forces to cancel themselves to some extent as N uj b is brought 
to 0. Letting f denote the resultant thrust force per unit mass of B , and letting fR, n 2 , and n 3 be three 
mutually perpendicular unit vectors fixed in N , one can obtain an accurate assessment of the issue by 
calculating three velocity increments as AV r = f f • n r dt (r = 1, 2, 3). In the present example, the total 
velocity increment is found to be 0.0073 m/s; using an initial mass of mo = (1.118 x 10 6 — 108) kg together 
with the rocket equation, one can determine the corresponding mass of SEP propellant to be Am = 0.2785 
kg. Thus, in this case, the failure to apply thrust as a couple is hardly detrimental. More generally, an 
upper bound on Am can be estimated by assuming f is constant in N for, say, 1800 s (30 min), with 
f = [(400hi + 400n2)/1.118 x 10 6 ] m/s 2 . Using the resultant of f, one obtains AV = -\/2(3.5778 x 10 -4 
m/s 2 ) (1800 s) = 0.91 m/s, and a corresponding mass of SEP propellant Am = 34.6 kg. In the worst case, 
then, the net savings in propellant mass is still almost 150 kg. 


7. Conclusions 

Changing the trajectory and rotational motion of an asteroid, whether for the purpose of retrieving it, 
deflecting it, or using it to deflect a larger asteroid, will likely involve some or all of the spacecraft motions 
considered here: approaching, hovering, synchronizing the rotational motion of the spacecraft with the 
asteroid, descending, ascending, and arresting the angular velocity of the asteroid. A preliminary study of 
the dynamics and control of these motions has been conducted in order to understand the magnitude of the 
required velocity increments, attitude control torque, and the associated propellant mass for a particular 
spacecraft and reaction control system. 

Gravitational forces exerted by the asteroid on the spacecraft are neglected. (Although this is justified in 
the case of a 10-m asteroid, it becomes less justified as asteroid mass increases.) Initial approach, hovering, 
descent, and ascent are examined for asteroids performing what is essentially simple rotational motion at 
angular speeds of 0.5, 1, 5, 10, and 60 revolutions per hour. The axis of rotation is essentially in the 
plane of the asteroid’s heliocentric orbit, and perpendicular to local vertical. Initial approaches that end in 
conditions suitable for hovering appear to be relatively inexpensive when using minimum-AV, multi-impulse 
trajectories. Closed-form expressions are presented for AV required to hover for one revolution. Hovering 
over one of the poles costs nothing; however, the cost escalates as the latitude approaches the equator, and 
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as the angular speed increases, and can become prohibitive. (Gravitational attraction may reduce the cost 
somewhat in the case of a large, slowly rotating asteroid.) As for descent and ascent, expense again increases 
in connection with equatorial regions and high angular speeds. Time of flight also enters the picture, and 
in some cases cost can be minimized with respect to time of flight. 

Two attitude control algorithms have been developed, and both are demonstrated with examples in- 
volving tumbling rotational motion of an unsymmetric rigid body. In the first case, the controller is used 
to synchronize angular velocity of the spacecraft with that of the asteroid, while the spacecraft’s axis of 
symmetry is kept parallel to a line fixed in the asteroid. The amount of propellant required is moderate, 
despite the asteroid’s high angular speeds. In the second case, an asymptotically stable controller arrests 
the angular velocity of an asteroid that is rigidly attached to the spacecraft. Rotational kinetic energy can 
be dissipated slowly or quickly, depending on the premium one places on control effort. A novel approach is 
applied to an example involving what is at first simple rotation: a perturbation induces tumbling, and this 
allows the use of long thrust moment arms that otherwise would not come into play. Propellant consumption 
is reduced by a factor of three, compared to an arrest in which simple rotational motion is maintained. 

In future work, we intend to improve the aforementioned analyses by accounting for gravitational attrac- 
tion of the asteroid. Solar arrays needed for the solar electric propulsion system will be large, and therefore 
solar radiation pressure will also be modeled. Treatment of the reaction control system will be refined by 
considering minimum firing times and minimum thrust levels. 
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